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Abstract 



We present a lattice QCD calculation oi B ^ ttIv semileptonic decay form 
, factors in the small pion recoil momentum region. The calculation is per- 

formed on a quenched 16^ x 48 lattice at /3 = 5.9 with the NRQCD action 
including the full 1/M terms. The form factors fi{y ■ k.,^) and f2{v ■ A;^) de- 
■ fined in the heavy quark effective theory for which the heavy quark scaling 

is manifest are adpoted, and we find that the 1/M correction to the scaling 
^ I is small for the B meson. The dependence of form factors on the light quark 

5h ' mass and on the recoil energy is found to be mild, and we use a global fit of 

the form factors at various quark masses and recoil energies to obtain model 
independent results for the physical differential decay rate. We find that the 
B* pole contribution dominates the form factor f^{(f') for small pion recoil 
energy, and obtain the differential decay rate integrated over the kinematic 
region > i8 GeV^ to be x (1.18 ± 0.37 ± 0.08 ± 0.31) psec-\ where 

the first error is statistical, the second is that from perturbative calculation, 
and the third is the systematic error from finite lattice spacing and the chiral 
extrapolation. We also discuss the systematic errors in the soft pion limit for 
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f^i^max) the present simulation. 
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I. INTRODUCTION 



The exclusive decay modes n^l^ui and —>■ p^l^v may provide us with the best 

experimental input to determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\Vub\- At present these decays are measured by CLEO with error of order 20%. A 
prerequisite for the determination of \Vub\ is an accurate calculation of the form factors 
involved in these semileptonic decays, but the theoretical prediction of the form factors for 
the entire kinematical range is still difficult. However, with the advent of the B factories, 
BaBar, Belle and CLEO III, we expect that the differential decay rate will be measured 
precisely as a function of momentum transfer in near future. This means that to determine 
\Vub\ we do not necessarily need the form factor for the entire kinematic region of g^, but 
calculations in a certain limited range of will practically suffice. 

Lattice QCD provides a promising framework to compute the form factors without re- 
sorting to specific phenomenological models. Exploratory studies have already been made 
by a few groups but more extensive studies are clearly needed to provide realistic 

predictions. In this work we attempt to compute the form factors and differential decay 
rates of 5 — > 'kIu for the momentum range q^ > 18 GeV^, which is set by the condition that 
the spatial momenta of the initial and final hadrons be much smaller than the lattice cutoff 
1/a, |k| -C l/a~ 2 GeV/c, to avoid discretization error. 

An important point in the calculation of the B meson matrix elements is to reduce the 
systematic error arising from a heavy quark mass M which is larger than 1/a. One approach 
adopted in the literature is to calculate the matrix elements with a relativistic action for 
heavy quarks around the charm quark mass and to extrapolate them to the bottom quark 
mass. Although this approach seems to work reasonably well in the recent studies of -B ^ nlu 
form factors @J^, the systematic error is magnified in the extrapolation and the heavy quark 
mass dependence would not be correctly predicted. This problem can be avoided by using 
a variant of the Heavy Quark Effective Theory (HQET), in which the the heavy quark is 
treated non-relativistically. 

A natural implementation of the idea of HQET on the lattice is the non-relativistic QCD 
(NRQCD) [^], which we employ in this work. With the NRQCD action the heavy quark 
mass dependence of the form factors can be reliably calculated , since the action is written 
as an expansion in terms of inverse heavy quark mass and higher order terms can optionally 
be included to achieve desired accuracy. In the B — >■ irlu decay near zero recoil of the pion, 
we find that the heavy quark expansion converges well at the next-to- leading order in 1/M. 

An alternative implementation of the HQET is the Fermilab formalism [^], in which 
results from the conventional relativistic lattice action are reinterpreted in terms of a non- 
relativistic effective Hamiltonian. This formalism shares an advantage similar to that of 



NRQCD, and has recently been applied to a 5 — nlu decay calculation |Tl 



In the application of the HQET to the B nlu decay, it is more natural to work with 
the form factors fi{v ■ A;^) and /2(f ■ k^^) |]12[, where v'^ is a heavy quark velocity and /c^ is 
a four-momentum of pion, rather than the conventional f^{q^) and f^{q^)- This is because 
the argument v ■ kj^, which is the energy of the pion in B meson rest frame, is well-defined 
in the limit of infinitely heavy quark mass, and the heavy quark scaling, i.e. fi.2{v ■ k-^) — *• 
constant as M — > cxd, is manifest in the new set of the form factors. 

We calculate /i,2(^ ■ ^tt) using the NRQCD action on a quenched lattice of size 16^ x 48 at 
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(3 = 5.9 corresponding to 1/a ~ 1.6 GeV. The action we use includes the full terms of order 
1/M. The 0(a)-improved Wilson fermion action is used for the light quark. We prepare a 
large statistical sample, accumulating 2,150 gauge configurations to reduce statistical noise 
which becomes large for states with finite momenta. This enables us to obtain good signals 
for the form factors for a finite spatial momentum of the pion. 

This paper is organized as follows. In the next section we briefly review the definition 
of the HQET motivated form factors fi,2{v ■ /Ctt) of Burdman et al. and their relation 
to the conventional form factors. We summarize the definition of the NRQCD action in 
Section and discuss matching of the heavy-light vector current on the lattice with that 
in the continuum in Section 0. We describe our lattice calculation in Section [V|, and 
the results are presented in Section Section [VI I| is given to a comparison with other 
lattice calculations, and phenomenological implications are discussed in Section |VII1| . Our 
conclusions are presented in Section |^ 



II. THE HQET FORM FACTORS FOR B ttLv 



The matrix element (7r(/c7r)|5'7^&|-B(p_B)) for the heavy-to-light semileptonic decay B 
ttIu is usually parameterized as 



{n{K)\qrb\B{pB)) = f^{q' 



+ nq')'^^^q'. (2.1) 



with pb and k^^ the momenta of the initial and final pseudo-scalar mesons and q 
When lepton mass is negligible, the momentum transfer q^ ranges from to 
m-j^Y. From the kinematics 



V ■ 



2 2 2 

— — q 



2m 



B 



Pb - k^. 

= (rriB - 



(2.2) 



where v = Pb/^b is the four- velocity of the initial B meson, a low q^ corresponds to a 
large recoil momentum of pion, for which the lattice calculation is not easy. In the other 
limit q^ ~ Q^axi however, the energy of the pion E^^ in the B meson rest frame is minimum, 
so that spatial momenta of the initial and final hadrons are small compared to the lattice 
cutoff, and the lattice calculation will give a reliable answer. 

In HQET, it is more natural to use and /c^ as independent four-vectors rather than 
and fc^. Burdman et al. defined the form factors fi{v ■ k^^) and f2{v ■ k^^) by 



{7r{k^)\qYb\Biv)) 



h{v-k^)v^' + h{v-k^) 



(2.3) 



where the heavy meson field is normalized with the factor 2f° instead of the usual 2p^, so 
that ^mB\B{v)) = \B{pb)). The new form factors are functions ofv-k-,^ and defined over the 
range [m^, (m^ — m^)/2m^]. As seen from definition ( [^.3D there is no explicit dependence 
on the heavy meson mass. Therefore, heavy quark scaling as M — oo is manifest, namely, 
fi,2{v ■ fcjr) become independent of M up to logarithms arising from the renormalization of 
the heavy-light current. Finite M corrections are given as a power series in 1/M. 
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The relation between the two definitions of form factors is given by 



f2{v ■ K) ^ fijv ■ K) 



V ■ K 



TUB 



(2.4) 



/ 



m 



B 



V ■ k. 



[h{v-K) + f2{v-K)] 



rriB 



fi{v ■ K) + 



ml 



{v ■ K) 



-J2{v ■ K) 



This indicates that f^{q^) and f^{q^) scale in the heavy quark limit as 

if f ■ /ctt is kept fixed. 

In the soft pion limit fc^ — > and ^ 0, we obtain simpler relations 



(2.5) 



(2.6) 
(2.7) 



f2{v ■ K) 

[fl{v-k^) + f2{v-k^)] 



(2.8) 
(2.9) 



from ( |2.4| ) and ( p.Sj) . The soft pion theorem implies that the scalar form factor f {q ) and 
the B meson leptonic decay constant fs are related as f^{qmax) = Ib/ fn, which means 



/l(0)+/2(0) 



fBy/m^ 



(2.10) 



The vector form factor f^{q^) may be evaluated using the heavy meson chiral lagrangian 
approach |T^, in which the B* pole contributes through a B*Btt coupling. One obtains the 
relation 



r r , ^ fB*^/rf}^ V ■ k^ 
Jm^/,(. . k.) = g + 



(2.11) 



where the vector meson decay constant fs* is defined by {Q\V^\B*{p)) = ifB*'m'B*^'^{p), and 
g denotes the B*Btt coupling. The B* propagator gives a factor l/{v ■ k^ + A^), in which 
= mB* — mB- Since the hyperfine splitting A^ ^ 46 MeV is much smaller than the 
'pion' mass, we consider in the lattice simulation that ( p^.llD depends little on on f ■ fc^. This 
behavior of /2 is actually found in our simulation. Equation (|2.11|) leads to the well-known 
vector meson dominance form for the form factor f^{q^) 



lim /+(g^ 



fB* 
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which is also reproduced in our calculation. 



2 ' 



(2.12) 
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III. LATTICE NRQCD 



We use the NRQCD formalism defined on the lattice to treat the heavy b quark 
without large discretization errors increasing as a power of aM. NRQCD is designed to 
approximate non-relativistic motion of heavy quark inside hadrons, and is expressed as a 
systematic expansion in some small parameter depending on the hadron considered. For 
the heavy-light meson system such as the B meson, the expansion parameter is given by 
^qcd/M, with Aqcd the typical momentum scale of QCD ~ 300-500 MeV. At the next-to- 
leading order in Aqcd/M, the lagrangian in the continuum Euclidean space-time is written 
as 



£cont ^ 
NRQCD — V 



Do 



2M 



9- 



(T B 

2M 



Q, 



(3.1) 



for heavy quark field Q represented by a two-component non-relativistic spinor. The deriva- 
tives Dq and D are temporal and spatial covariant derivatives respectively. The leading order 
term Dq represents a heavy quark as a static color source. The leading correction of order 
^qcd/M comes from Z)^/2M, which gives a non-relativistic kinetic term of heavy quark. 
Another contribution of order Kqcd/^ is the spin-(chromo)magnetic interaction cr- B/2M, 
where B denotes a chromomagnetic field strength. In the usual HQET approach, only the 
leading terms are present in the effective lagrangian and corrections of order Aqcd/M is 
incorporated when one evaluates a matrix element (O) of some operator O by including 
terms such as {TO J d^xQ^{D'^ /2M)Q). In contrast, in the NRQCD approach we include 
the correction terms in the lagrangian (|3.1|) and evaluate the matrix elements with the heavy 
quark propagator including the effect of order KqcdI^- 

An important limitation of the NRQCD lagrangian is that the heavy quark expan- 
sion is made in the rest frame of a heavy quark. Since the expansion parameter is p/M, 
where p is a typical spatial momentum of heavy quark, the lagrangian is valid only in the 
region where the heavy quark does not have momentum greater than 0{Kqcd)- Therefore, 
in the study of the heavy-to-light decay, the momentum of initial B meson must be small 
enough. Although it is possible to construct the action expanded around a finite heavy quark 
velocity, the heavy quark velocity is renormalized by radiative correction since the lattice 
violates Lorentz symmetry |T^,|T5|, which gives rise to an additional important systematic 
corrections. We, therefore, do not take this strategy and consider the discretization of the 
Lagrangian ( |3.1| ). 

The lattice NRQCD action we use in this work is 

Snrqcd = - KQ{x,y))Q{y) + ^x\x){6^^y - K^{x.y))x{y). (3.2) 



In addition to the non-relativistic heavy quark field Q, we write the term for the anti-particle 
field X for completeness. The kernels to describe the time evolution of heavy quark are given 
by 



Kqix^y) 
Kx{x,y) 



1 



1 - 



~2n 



aSH' 



si'^ul 1 



2n 



5r'U, 1 




(x,?/), (3.3) 
x,y), (3.4) 
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where n denotes a stabilization parameter introduced in order to remove an instability 
arising from unphysical momentum modes in the evolution equation The operator ^4^^ 
is defined as 6'^\x,y) = Sx^±i,y4^Sx,y, and Hq and 6H are lattice Hamiltonians defined by 

where A*-^-* = X]?=i '^i^'* is a Laplacian defined on the lattice through a\^\ the second 
symmetric covariant differentiation operator in the spatial direction i. In (|3.6| ) the chromo- 

magnetic field B is the usual clover-leaf type lattice field strength M. In these definitions, 

(2) 

the lattice operators A^^ and B are dimensionless, i.e. appropriate powers of a are under- 
stood. The space-time indices x and y are implicit in these expressions. The bare heavy 
quark mass Mq is distinguished from the renormalized one M. 



The lattice action (|3.2| ) describes continuum NRQCD ( p.l|) in the limit of vanishing lattice 
spacing a at tree level. In the presence of radiative correction, however, power divergence 
of form a"/(aMo)™' with positive integers n,m can appear. This is due to the fact that 
NRQCD is not renormalizable, and the action should be considered as an effective theory 
valid for small l/(aMo). This means that the parameters in the lattice action ( |3.2| ) should 
be tuned to reproduce the same low energy amplitude as the continuum QCD up to some 
higher order corrections. One may use perturbation theory to achieve this tuning. For 
example, one-loop calculation of energy shift and mass renormalization was carried out for 
lattice NRQCD by Davies and Thacker and by Morningstar [|T^ sometime ago, and then 
by ourselves [p^^ -pOj for the above particular form of the NRQCD action.[] To improve the 
perturbative expansion we utilize the tadpole improvement procedure where all the gauge 
links in the action ( p.2| ) are divided by its mean field value Uq determined from the plaquette 
expectation value as Uq = ((TrUp)/3)^/^. This tadpole improvement will give rise to 0{g^) 
counter terms in the Feynman rules. The one- loop tuning of the coupling constant in 
front of the spin-(chromo)magnetic interaction term (|3.6| ) has not yet been performed. We, 
therefore, use the tree level value cs=l after making the tadpole improvement. 

The relativistic four-component Dirac spinor field h is related to the two-component non- 
relativistic field Q and x appearing in the NRQCD action ( |3.2| ) via the Foldy-Wouthuysen- 
Tani (FWT) transformation 





'-11-5^1 I I. (3.7) 



where V is a symmetric covariant differentiation operator in a spatial direction. 



^ We note that the evolution kernels ( |3.3D and (|3.4| ) are slightly different from the definition used, 
for example, in |17|, where the (1 — ai/o/2n)" terms appear inside of the (1 — a6H/2) terms. 
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IV. MATCHING OF THE HEAVY-LIGHT CURRENT 



Since we use the lattice NRQCD action of the previous section, the continuum heavy- 
hght vector current 97^6 in ( |2.1| ) must be written in terms of the corresponding operator 
constructed with the lattice NRQCD heavy quark field h. This matching of the continuum 
and lattice operators has been calculated using the one-loop perturbation theory by Morn- 



ingstar and Shigemitsu pl|j2^ . In this section we summarize their results and specify our 
notations. 

In the one-loop matching of the continuum operator to the lattice operators, we have 
to consider dimension-four operators in addition to the leading dimension three operator 
57^/;., in order to remove the error of order as^^Qco/ M and asaKqcD- The former is the 
radiative correction to the FWT transformation (|3.7|) and the latter appears in the 0(a)- 
improvement of the lattice discretized operator. Thus the following operators are involved 
in the calculation. 



(0 



(1 



(2 



(0 



(1 



(2 



(3 



(4 



974/1, 



2aMo 
1 

~2aMo 
1 

2aMo 
1 

2aMo 
1 



9747 ■ V/i, 
97- V h, 

q'jkl ■ V/i, 
97- V jAlkh, 



2aMQ 



2aM, 



(4.1 
(4.2 

(4.3 

(4.4 
(4.5 

(4.6 

(4.7 

(4.8 



The heavy quark field h is obtained from the two-component field Q through the FWT 
transformation (|3.71 ) .p| For the light quark q we employ the 0(a)-improved Wilson fermion 



The one-loop matching is given by 
1 



•f/cont 
^4 



1 + as 



.(0) 
14 



cont 



-ln(aMo)+pi;^ 

TT 

1 - ' - - (0) \ T.(0) 



(i)t/(i) 

k 



+Ois 



-ln(aMo)+pg 

TT 



(3)t.(3) 



-i^ln(aMo)+pg 



(4) 



(4.9) 



(4.10) 



In the definition used in |22| the heavy quark field before the FWT transformation [Q 0)^ 



appears in the definition of operators. Matching coefficients for 
when we use the above definition. 



(1) 



and Vi^^ must be converted 
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and the numerical coefficients py] and py^ are summarized in Tables | and || for several 
values of aM^. 

As we mentioned earlier, the NRQCD action employed in this work is slightly different 
from that of Morningstar and Shigemitsu p2|]. We have therefore independently calculated 



the wave function renormalization and the vertex correction for the temporal component V4, 
and found that the difference of the finite constants p's between the two actions is small, 
e.g., ~ 4-9% for the vertex correction. Therefore, for the spatial vector current, for which 

assuming 



the one-loop calculation with our action is missing, we adopt the coefficients of |]22[ 
that the error is negligible. In Table | the results of our calculation for py] are listed, while 



results for py^ in are interpolated in aMo and given in Table || for our parameter values. 



V. LATTICE CALCULATION 
A. Lattice setup 

Our quenched lattice calculation is carried out on a 16^x48 lattice aX {3 = 5.9 with the 
standard plaquette action for gluons. The inverse lattice spacing 1/a determined from the 
string tension is = 1.64 GeV. The scaling violation has been found to be small for our 
choice of the heavy and light quark actions over l/a~ 1 — 2.5 GeV in the heavy-light decay 
constant |T^. 

The parameters we choose for the heavy and light quarks are a subset of those simulated 
1^. We take four values of the bare mass aMo, 1.3, 2.1, 3.0 and 5.0 for the heavy quark. 



m 



over a range of the physical heavy quark mass between 2 and 8 GeV. The stabilization 
parameter n is set to 3 (for aMo=1.3 and 2.1) or 2 (for aAfo=3.0 and 5.0) so as to satisfy 
the stability condition n > 3/ (aMo). We use the 0(a)— improved Wilson action for the light 
quark with the clover coefficient Csw=l-580, which is evaluated at one-loop with the tadpole 
improvement. Four values 0.13630, 0.13711, 0.13769 and 0.13816 are chosen for the hopping 
parameters in our simulation, where the critical value Kc is 0.13901. 

We accumulate 2150 quenched configurations to reduce the statistical error for matrix 
elements with finite spatial momenta. Each configuration is separated by 1000 pseudo-heat- 
bath sweeps after 10 000 sweeps for thermalization and fixed to the Coulomb gauge. As we 
will see, even with this large number of statistics, signals for heaviest heavy quark or lightest 
light quark are not clean enough to extract the ground state. 

B. Correlators 

The form factors are extracted from measurements of three-point correlators 

C-'^^''''\u,t,tB;q,PB) = E e-^1-''e^PB-y{O'Au,0)Vi^^\t,x)O'J{ts,y)), (5.1) 

x,y 



of the vector currents (4^)^(4^) with the initial B meson and daughter pion interpolating 
fields and O^, respectively. The interpolating fields are defined by 
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0^{t,x) =q{t,x)-f5q{t,x) 



J2M\rMt,x + r) 



75 



J2M\r'\)q{t,x + r') 
r' 



0^it,x) = q{t,x)-f5h{t,x) 



Ol{t,x) = g(t,x)75 



J2M\rmt,x + r) 



(5.2) 

(5.3) 

(5.4) 
(5.5) 



where the operators with superscript L represent a local field, while the smeared operators 
defined on the Coulomb gauge fixed configurations are labeled with S. The smearing func- 
tions (f)i{r) and </'h(r) are parameterized by (f)i{r) = exp(— a^r''') and 0/i(r) = exp(— a/^r^''), 
with parameters ai h and bi h determined from a measurement of light-light and heavy-light 
meson wave functions |19|. The wave function of the light-light meson (j)i{r) is almost in- 
dependent of the light quark mass, and we use {ai,bi) = (0.27,1.13). The wave function 
describing the spread of heavy-light meson (f)h{f ), on the other hand, depends significantly 
on the heavy quark mass, i.e. becomes larger as heavy quark mass increases. The 



numerical values of {a^, b^j are given in Table pi . 



The B meson interpolating field is fixed at the time slice ^^=24. The light quark 
propagator corresponding to the spectator quark is solved for a smeared source at t^, and the 
source method is used at the time shce t,r=0 to obtain the daughter light quark propagator 
with momentum insertion —q. The heavy-light current is then constructed at t, which is in a 
region t^, < t < tsi with the daughter light quark propagating from and a heavy anti-quark 
evolving back from Ib- Another momentum is inserted at t. With this combination of 
momenta, the initial B meson has momentum Pb and the final pion travels with momentum 
fcvr = Pb ~ since the fixed source at emits a heavy-light meson with any momentum. 



The momentum combinations measured in our simulation is summarized in Table IV. Since 



the statistical noise grows exponentially as exp((i?(p^) — E{0))t) for the finite momentum 
state with energy E{p'^), the spatial momentum one can measure with reasonable signal is 
rather limited. In fact, even in our high statistics data, the maximum momentum we could 
take is (1, 0, 0) in unit of 27r/La as we shall discuss in the following sections. 

The three-point function (|5.1| ) is dominated by the ground state contribution for large 
enough separation of operators <^ t <^ Ib 



C^'''^''^\u,t,tB;q,PB) 



ZUpb) 



2E^{K)2Eb{pb) 
X exp [—ET,{kT,){t 



{AKW^^WipB)), 

- t^) - Eund{PB)itB - t)) . 



(5.6) 



The overlap amplitudes Z'^(fc^) and Z^(pb) of the interpolating operators with the corre- 
sponding ground state are evaluated from the two-point correlators defined by 

C-'-\u,t;K) = 5:e'^-^(Of(t.,0)Of (t,^)) ^ giM!e-^^(fc^)(*-*^), (5.7) 

C^'^^t^tB-.PB) = Y.^-'P--''{Ol{t,x)Ol\tBM ^ M(P^e-^-(P-)(*-*), (5.9) 
X -^^bkPb) 
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C^'^'^t^tB^PB) = Y.^-"^--''{O^BMOl\tBM ^ ^%l^%le-^--(P-)(*-"*). (5.10) 
X ^Eb[Pb) 

The ground state energy of the heavy-hght meson -EbindlPB) represents a "binding energy", 
as the bare heavy quark mass is subtracted in the NRQCD formahsm. In the state normal- 
ization in ( |5.6| ) and in ( |5.9| )-( |5rT0| ), on the other hand, the heavy-hght meson energy Eb{Pb) 
including the bare heavy quark mass enters in the denominator. 

In practice, we calculate a ratio i?^** of the three-point and the two-point 

functions 



C-'-\K.t-K)C^'^^'{t,tB-.p^) ^ ZL{k^)Zj^{ps) 



which becomes constant in the asymptotic limit. The overlap amplitudes with the smeared 
interpolating fields Z^lk^^) and Z^{p^) cancel between the numerator and the denominator. 

Typical examples of the ratio R^i^ {t]k.,^,p^) is plotted in Figure |l], in which the data at 
K=0. 13711 and aMo=3.0 are shown for five choices of the momentum combinations. For all 
these plots we find clear plateau in the large t region, where the current is closer to the B 
meson interpolating field than to the pion. This is due to the fact that the smearing with 
the measured wave function works better for heavy-light than for light-light meson, because 
the heavy-light system is less relativistic and the description with the wave function is more 
appropriate. The fit result is indicated by horizontal lines. 

The data become noisier for lighter light quark masses with a fixed heavy quark mass, 
or for heavier heavy quark mass with fixed light quark mass. As a result, we are not able 
to extract signals for our lightest light quark k=0. 13816, except for a few cases when the 
daughter pion does not have finite spatial momentum. We also note that we carried out 
simulations for one additional heavy quark mass aMo=10.0. We found, however, that the 
signal is intolerably noisy, so that we do not use those data in our analysis. 



C. Matrix elements 



In order to obtain the matrix element (7r(/c^)|V^^*''|i?(pB)) /^from (|5.11| ), we have to elim- 



inate in the denominator. For this purpose we fit the smeared-smeared and 

smeared-local two-point functions with a single exponential as in (|5.71) and (|5.8|) for the ex- 



traction of Z^{k^)/^E^{k^), and (|0| ) and ( ^.IQ ) for Z^{p^)/ ^ Eb{Pb)- We then obtain 
a combination 

V.<.(...p,).<«l£M, ,5,12) 

^E^ik^)EBiPB) 

Numerical results are listed in Tables |V|- |V111| for each light and heavy quark mass. The first 
column denotes the momentum configuration as shown in Table |V[ 



VI. RESULTS FOR THE FORM FACTORS 
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A. Energy-momentum dispersion relations 



In order to extract the form factors from the matrix elements ( [5.12D , we have to determine 



the meson energy of initial and final states for given spatial momenta. It may be obtained 
either by assuming a continuum dispersion relation or by actually measuring the meson 
energy with the given momenta. 

For the pion, which is relativistic, the continuum dispersion relation is written as 

E^{K)^ = M^ + kl. (6.1) 

The measured values of {aET,{kT,)y for momenta fc^ = (1,0,0) and (1,1,0), in unit of 
27r/La, are given in Table ^ and also plotted in Figure ^ for each light quark mass we 
calculated. We find a nice agreement with the expectation ( |6.1|) . The relation ( |6.1|) may 
be modified on the lattice due to lattice artifacts; a possible form is given by replacing akj: 
with sin(afc^), which satisfies the periodic boundary condition. The magnitude of such an 
effect is not significant, though, since the momentum considered is small enough and the 
difference between ak^, and sin(aA;^) is less than 3%. 

The dispersion relation for the heavy-light meson is well described by the non-relativistic 
form 

EUPb) = ^bin(O) + (6.2) 

in which the meson mass Mb appears in the kinetic energy term.0 In NRQCD, the heavy- 
light meson mass is written in terms of the bare mass gMq and the binding energy aEun{0) 

as 

aMe = Z^aMo - aEo + aEun{0), (6.3) 
where qEq is an energy shift and is a mass renormalization factor. Both factors are 



calculated at the one-loop level [|16],|T^|T9 



aEo = asA, (6.4) 
Zm = l + asB, (6.5) 

and the numerical coefficients A and B are given in Table I of |19[. The heavy-light meson 
mass evaluated with ( |6.3|) using the V^-scheme couphng av{q*) p4| at g* = 1/a is listed 
in Table Since the one-loop correction partially cancels between Z^aMo and aEo, the 
uncertainty due to the choice of q* is small, i.e. at most 3% for aMo=1.3 and even smaller 
for larger aMg. 

In Figure |^, a comparison is made of our simulation data with the form of Eq. ( |6.2|) 
in which the value of Mb evaluated according to Eq. ( |6.3| ) is substituted. We find a good 
agreement except for the data at fi;=0. 13630. Even in the worst case, the disagreement does 



^ Here we use a capital symbol Mb to represent the generic heavy-light meson mass we deal with 
on the lattice, whereas keeping tub to denote the physical B meson mass. 
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not exceed 1%. Therefore, we employ the dispersion relation (|6.2| ) with the perturbatively 
estimated meson mass qMb in the following analysis of the form factors, rather than using 
the measured binding energy, which has significant statistical errors and complicates our 
analysis. The same strategy is taken for the pion energy, namely we use the relation ( |6.1|) 
with the measured value for aikL. 



B. Form factor extraction 



The continuum matrix element is obtained from V^^^^H^tjPb) defined in (|5.12| ) using the 
matching formula of the vector current ( [4.9| )-( ^?T0D as 



^4 (fc.,PB)= l+« 



■ln(aMo)+pS^^ 



TT 



(i)t>(i), 

4 



(6.6) 



(^tt; Pb) 



1 + as 



(0) 



-HaMo)+p'" 



TT 



vt'\K,PB) + aJ^!vt'\K,PB) 



+Ols 



-a. 



- In(aMo) + pg 



vt'\K^PB) + c^sP'v!vr{K,PB 



(3)t>{3), 



-— ln(aMo)+pl?^ 



vt\k.,PB). 



(6.7) 



We use the F-scheme coupling av{q*) for the coupling constant as- Since the scale q* which 
dominates the lattice one-loop integral is not yet known, we examine the uncertainty in the 
scale setting by calculating the form factors at g* = 1/a and at vr/a. We use the difference in 
the results, which is the two-loop effect of O(a^), as an estimate of higher order perturbative 
errors. The numerical value of the coupling is av'(l/a)=0.270 and a;y(7r/a)=0.164 at /?=5.9 
in the quenched approximation. 



From the definition of fi{v ■ fc^) and f2{v ■ fc^) given in (2^), we obtain the following 
formula for the form factors 



IE^{K)Eb{pb) T rcont 



AM, 



B 



f2{v ■ K) 



1 - 



Ml 



(v ■ K) 



I E^{K)Eb{P b). 



vr\K,PB) 



(6.8) 



, (6.9) 



where = {Eb{pb),Pb)/^b and fc^ = (-E'^(A;^), fc^). By construction, for the initial B 
meson at rest, /i(f ■ A^tt) + /2('J^- ^tt) is proportional to the temporal component V2^°"^*(fc7r, Pb); 
while /2(f ■ /Ctt) comes from the spatial component I4™"*(fc7r, Ps). Even for a B meson with 
momentum (1,0,0), the velocity is small {pb/Mb — 0.07-0.2 depending on the heavy quark 
mass), and the major effect is from the temporal or spatial component for fi{v ■ kj^) + f^iv ■ k^^) 
or for /2(f ■ /ctt), respectively. 

An example of the form factors is plotted in Figure ^ for aMo=3.0, which is close to 
the b quark mass, and fi:=0. 13630. The point of smallest av ■ k^^ corresponds to the zero 
recoil configuration, i.e. the initial and final particles are at rest so that av ■ k^^ = aM^^. At 
that point, only the temporal component V2^°"*(fc7r, Pb) can be measured while the spatial 
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component vanishes. The momentum configuration p^=(l,0,0) and fc7r=(0,0,0) gives a very 
similar av ■ /c^, because of large heavy quark mass and small spatial velocity. As a result, the 
data point almost lies on top of that at zero recoil. We are not able to measure f2{v ■ /c^) 
reliably at this point, since the value of the spatial component V^°^^{k.„,p^) is too small. 
There are four other momentum configurations (see Table [TVD , for which both fi{v ■ kj^) + 
f2{v ■ A;^) and f2{v ■ kj^) are measured. Among them, two momentum configurations sharing 
the same fc^=(l,0,0) and having different have almost identical values of av ■ /c^ for the 
same reason as above, and cannot be distinguished from each other in the plot (the middle 
point of the three filled data points). 

From Figure ^ we also see that the effect from choosing Q;y(l/a) (circles) or ay [it /a) 
(squares) is small; it is smaller than the statistical error except for the zero recoil point 
where the statistical error is minimum. Therefore in the following analysis we use the data 
with ay(l/a). In the final results we will include their difference in the systematic error 
estimation. 



C. Heavy quark mass dependence 

As we discussed in Section |T|, the heavy quark scaling is manifest for the form factors 
/i(f ■ k^,) and /2(f ■ A;^). Namely, fi^2{v ■ k^,) behaves as a constant at the leading order of the 
1/M expansion. Here we examine the heavy quark mass dependence of fi{v ■ /c^) + /2(w ■ k-,^) 
and /2(f ■ /Cjr) explicitly by comparing the results with different heavy quark masses. 

In order to remove the logarithmic dependence on the heavy quark mass that appears 
from the matching of the vector current between the full QCD and lattice NRQCD ( |4.9| )- 
( [4.10| ), we define the renormalization group invariant form factors $1+2 (^ ■ k.,^) and (^2{v ■ k^^) 
as 

2//3o 

[fi{v ■ k^) + f2{v ■ k^)] , (6.10) 

. . , ..v^ (6.11) 

where (3q denotes the first coefficient of the QCD beta function. We note that Mb is the 
heavy-light meson mass measured on the lattice for a given oMq, while is the physical 
B meson mass. 

Figure ^ shows $i+2(f ■ k^,) and ^2{v ■ k-,,) for several values of oMq. We find that the 
1/M correction gives only a small effect in the range of the heavy quark mass we explored 
which correspond to 2-8 GeV. In fact, there is no significant shift in the magnitude of the 
form factors by the change of the heavy quark mass. A small effect can be seen in the value 
of av ■ kjr for a couple of momentum configurations for which ■ k^^ ^ 0. However, it does 
not seem to change the global shape of the form factors. 

The small 1/M correction we found is of great phenomenological importance, as it jus- 
tifies the use of heavy quark symmetry to predict the B — > nlu form factors from that 




of -D ^ nlu and D Klv ||T2[. We discuss this method and possible uncertainties in 
Section IVTm. 
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D. Light quark mass dependence 



In order to obtain the physical form factors we need to extrapolate our result to the 
physical light {u and d) quark mass. For this purpose we examine the light quark mass 
dependence of the form factors fi{v-kT,) + f2{y-k.„) and /2(f ■A;,^) using the data, which covers 
a range 0.45-0.80 GeV of the pion mass. Unfortunately, the signal is badly contaminated 
by statistical noise for the lightest data, so that we are not able to extract the form factors 
except for the zero recoil limit of fi{v ■ fc^) + f2{v ■ fc^). For other three k values, the data is 
fully available and we mainly use them to see the light quark mass dependence. 

Figure ^ shows the measured form factors at four different light quark masses. The 
heavy quark mass is fixed at aMo=2.1. Since the minimum value of av ■ k.^ is aM^r, the range 
of av ■ kj^ where the data is available moves to the left hand side as the light quark mass 
decreases. On the other hand, change of the value of the form factors fi{v ■ k^,) + /2(f • k^,) 
and /2(f ■ /Cjr) is not significant if we compare the data for a given momentum configuration. 
For instance, the values of a^^'^[fi{v ■ k^) + f2{v ■ k^,)] stay almost constant around 0.68 over 
the range v ■ k-j^ = 0.27 - 0.49, which correspond to the lightest and the heaviest data. If we 
look at the change at fixed v-k^^, there is an apparent downward shift of /i (f ■ ^tt) + /2 ("^ ■ ^tt) ■ 
This is due to a negative slope in av ■ k^^ in the data at fixed light quark mass. On the other 
hand, for f2{v-k.„) the light quark mass dependence is less significant, since the data at fixed 
K does not seem to have a non-zero slope. 



E. Global fit 

In order to extract the physical form factors, we have to consider the dependence on three 
parameters, i.e. the inverse heavy meson mass 1/Mb, the light quark mass and the energy 
release v ■ k^^. The heavy quark effective theory together with the chiral perturbation theory 
suggests that we can expand the form factors in powers of 1/Mb and rrig. On the other hand, 
there is no theoretical guide in the functional dependence on v ■ k^^. Therefore, in fitting the 
data we use a Taylor expansion around an arbitrary chosen point v ■ kj^ = [v ■ k.j^)o, which in 
practice we take in the middle of the measured range. Thus we employ the following form 
to fit the data, 

^(100) 

aV2$i+2(^ ■ k^) = Cf^f + ^ + Cl'i^^am, + (cf!^^^ + Cf^^^am,) [av ■ k^ - {av ■ k^ 
+Cf^2 [o^ ■ K - {av ■ k^)Qf , (6.12) 

^(100) 

a^lH2{v ■ k^) = Cr°) + + Cf'^^am, + C^'^ [av ■ k^ - {av ■ k^)o] , (6.13) 



aM 



B 



where the superscript {ijk) for the coefficient denotes the order of expansion in l/aMs, arriq 
and [av ■ A;^ — {av ■ /c,r)o]; in the given order. The fit results for {av ■ /i;^)o=0.5 is listed in 
Table IXm 



The choice whether we keep a certain term in ( |6.12 )-( |6.13 ) or drop it is empirical. 



Our experience in the calculations of the heavy-light decay constant and the B parameters 
suggests that both the 1/Mb and arriq expansions can be safely truncated at the first order. 
This is consistent with an argument of naive power counting assuming that the relevant 
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mass scale is around ^^qcd- We find that that is indeed the case also for the B ttIu form 
factors as we shall discuss in the following. 

In ( |6.12D -( |6.13| ) the l/Mg expansion is truncated at first order, since the 1/Mb correction 



is not significant as discussed in Section [VI C| so that there is no sensitivity to higher order 



corrections. Even the first order coefficients C1+2''' ^ind C2^^^^ are consistent with zero within 
the statistical error. This truncation is also consistent with our NRQCD action, because we 
do not include the correction of order 1/Mq or higher in the action ( |3.2| ). 

The light quark mass dependence of a^/^$i+2(f ■ /Cjr) is consistent with a linear function 
if we fix av ■ k^^ at [av ■ ^7^)0=0. 5, for instance. Thus we truncate the expansion in arug at the 
first order. We also keep a cross term with the leading [av ■ fc^ — {av ■ fc7r)o] correction, but its 
coefficient C[^j^2^ is consistent with zero. For a^/^$2('^ ■ ^tt) the light quark mass dependence 
is not significant as discussed in Section | V 1 D|. Although we keep the first order correction 



to be conservative, its coefficient 6*2°^°'' is almost consistent with zero. 

As for the functional dependence of the form factors on av ■ k.„, we include the [av ■ k,, — 
(af ■/c^)o]^ term for a^/^$i+2(f -fcyr), while the second order term is neglected for a^^'^^2{v-k^). 
The reason is that we find a significant slope in a^/^$i+2('y ■ k^^), so that a higher order term 
[av ■ /ctt — (of ■ A;,r)o]^ is also included for safety. The other form factor a^^'^^2iv ■ kj^) behaves 
almost like constant, and it is enough to keep the first order term. 

While we introduce several terms for which the coefficient is not well determined, i.e. 
consistent with zero, this does not mean our results for the physical form factors have large 
uncertainty, as far as we use the results for an interpolation in the relevant parameters. 
For example, the heavy quark mass we simulate covers the b quark mass, so that only an 
interpolation is required. For the parameter [av ■ k.,, — {av ■ /c7r)o], we restrict ourselves to 
consider the region where the data is available. Therefore, we can obtain the physical form 
factors in the region 0.67 GeV < v ■ k^^ < 0.96 GeV reliably. Outside this region, the fit 
( |6.12| )-( |6.1^j| ) appears to introduce a large uncertainty. For the light quark mass, we have 



to consider an extrapolation to the physical limit of u and d quarks. This increases our 
statistical error significantly. 

The fit results are shown in Figure |^ (heavy quark mass dependence) and in Figure |^ 
(light quark mass dependence). In Figure ^ we also plot the limit of physical light quark 
mass (thick curves), which is obtained by setting aniq to the physical average up and down 
quark masses in (|6.12|) -( |6.13|) . 

The form factors f\{v ■ k-^) + /2(f ■ k-,,) and /2(f ■ A;^) for the physical B — > nlu decay is 
plotted in Figure |^. The region where the lattice data is interpolated in [av ■ /c^r — {av ■ fc7r)o] 
is plotted with symbols. Going outside of that region requires an extrapolation, so that the 
error shown by dashed curves rapidly grows. 



F. Soft pion theorem 

In the soft pion limit, i.e. and kj^ 0, the following relation ( 2.10| ) holds 



n<lL.) = -4^ [/i(0) + /2(0)] = ^. (6.14) 

It is an important consistency check to see whether one can reproduce this relation in the 
lattice calculation. 
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In Figure ^ we plot the result of the fit ( |6.12|) by an open triangle and compare it with 
the lattice calculation of fs/ fn (filled triangle) [^. The data is presented at fixed heavy 
quark mass aMo=3.0. We should note that the soft pion limit in ( |6.12| ) is far from the region 
where fi{v ■ k.^) + f2{v ■ k.^) is obtained by interpolation. Therefore, we expect substantial 
systematic uncertainty in the fit result. In fact, Figure |^ demonstrates that the extrapolation 
to f ■ = is not yet very stable. 

The soft pion limit can also be achieved along a fixed momentum configuration. Namely 
we may extrapolate the data for each light quark mass at zero recoil. In this case, however, 
the momentum transfer v ■ fc7r(= ^-n) changes during the extrapolation, so that we have to 
consider a fit with two terms amq and aM^^.^ Because of the PCAC relation oc mg, it 
means a quadratic fit in ^airiq. We plot two extrapolations in Figure ||: a linear form in 
amq (dashed line) and a quadratic fit in ^arnq (solid curve). Although the effect of the term 
^auiq seems very small in the data and is only seen at the lightest quark mass, it raises the 
soft pion limit for the quadratic fit. The result is consistent with the global fit ( |6.12 ). Thus 
we consider that the disagreement of f^iq^ax) with fs/ f-w, which seemed to be a serious 
problem by only looking at the naive linear extrapolation with only data at zero recoil, is 
in fact more of a problem in the subtle chiral extrapolation or in the model uncertainty of 
momentum extrapolation. 

The UKQCD [^,^] and APE |7,27] collaborations found in their study with the rela- 
tivistic heavy quark action that the soft pion relation ( |6.14| ) is satisfied. It should be noted, 
however, that their method of chiral extrapolation corresponds to our "global fit" method, 
and the measured kinematical region is far from the soft pion limit. Therefore the result in 
the soft pion limit should depend on how the extrapolation is made. They employed a pole- 



like model for their fit function. Thus their results in the soft pion limit contain some 
uncertainty which is not well controlled just as ours, although the results in the kinematical 
region obtained by interpolating the lattice data do not suffer from such uncertainties. 

Judging from the size of the uncertainties it is too early to consider the deviation from the 
soft pion relation as a serious problem. This problem can be studied with much statistically 
significant data with larger number of momentum points and light quark masses so that the 
extrapolation mv -k^^ towards the soft pion limit becomes more stable, which is still beyond 
the scope of this paper. 



G. Pole dominance 

In the soft pion limit, the heavy meson effective lagrangian predicts the B* pole domi- 
nance ( f^.llD , that is 



lim /,(. ■ k.) = r\ . (6.15) 



^ As discussed in , one should include a term which is linear in oMt^ when v • /ct^ (or in the 
relativistic form) varies during the extrapolation of form factors. The fit becomes more stable if 
one first interpolates to a fixed v ■ k-,^ (or q^), and then extrapolate in aruq. Our strategy to employ 
the global fit (|6.12| )-( |6.13| ) is equivalent to this method. 
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Since the hyperfine splitting = Mb* —Mb is much smaller than the momentum transfer v- 
kir we measure, we can approximate its functional form by a constant in our data region. Our 
data support the constant behavior and give g{fB^^/omB^/2f^,) = 0.35(18) which reduces 
to 5^ = 0.30(16). This agrees with the phenomenological value extracted from D* Dir 
decay 0.27(6) ||29|, and also with the recent lattice calculation g = 0.42(4) (8) ||3^, which 
is obtained for the static heavy quark. The agreement suggests that the 1/M correction is 
small for the form factors. 



H. Systematic Errors 

We now discuss possible sources of systematic errors and their estimates. Since the 
statistical error, the discretization error of O(a^), the perturbative error of O(a^), and the 
chiral extrapolation error are large, we only consider these dominant sources of errors and 
neglect other subleading errors such as 0(a^/(aM)), O^a'^aAgcD), O^asAQco/M) and so 
on. 

The size of the two-loop order correction is only known by explicit computation, which 
is beyond the scope of this paper. Instead, we estimate the size of the perturbative error of 
0(q;^) as half of the difference of values for q* = ir/a and 1/a. The typical sizes are 1.5% 
for /i(f ■ k^,) + f2{v ■ kj,) and 3.5% for /2(f ■ kj,). The reason for the error for /2(f ■ k^,) being 
larger is that the one-loop renormalization coefficient for heavy-light vector current is larger 
for the spatial component than for the temporal one and that the matrix element of the 
spatial component gives larger contributions to /2(f ■ /ctt) in the small recoil region. 

The discretization errors of ©(a^Ag^^) and of O(a^fc^) are also important. The former 
error is common to most lattice simulations using the 0(a)-improved actions, and through 
an order counting we estimate it to be 3% sX j3 = 5.9 assuming that the typical momentum 
scale Aqcd is around 300 MeV. The latter is specific to the present work since the error due 
to nonzero recoil momenta appears only in the study of form factors. As the pion momentum 
treated in our calculation is at most 27i/L (L=16) in the lattice unit, we estimate this error 
to be about 16% using the order estimation. 

The error in the chiral extrapolation is another major source of the systematic error. 
Since we have data at only three k values except for the zero recoil point, it is not practical 
to test different functional forms in for the chiral extrapolation. We instead estimate the 
corresponding error in the form factors by taking the square of the difference between the 
result of the chiral limit and that of the lightest k. This gives 10% for fi{v ■ kj^) + f2{v ■ k^^) 
and 1% for /2(f ■ k^^). 

The total error is estimated by adding these errors in quadrature together with the 
statistical error. In Figure |^ the form factors fi{v ■ k.„) + /2(f ■ A;^) and f^iy ■ fc^) are plotted 
with the estimated systematic uncertainties. Numerical results are listed in Table PCIII] . 



VII. COMPARISON WITH OTHER CALCULATIONS 
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A. fi{v ■ k^) and f2{v ■ K) 

The Fermilab group calculated the form factors at the h quark mass using the non- 
relativistic interpretation of the relativistic lattice action ||lT|. A comparison is made with 
our results for the HQET form factors fi{v ■ k^^) + f2{v ■ A;^) and f2{v ■ k.,^) at the same (3 



value employed, /?=5.9, in Figure |T0[ We find a reasonable agreement for f2{v ■ k^,), but for 
fi{v ■ fc^) + f2{v ■ kj,) our data seem substantially lower. 

Since both NRQCD and Fermilab action are the two variants of the nonrelativistic effec- 
tive action, there should be no fundamental difference in the result. There can be, however, 
two possible reasons for the disagreement. One is the difference in the renormalization fac- 
tor. The other is the difference in various systematic errors which arise from the choice of 
parameters such as the lattice size, smearing methods, fitting procedures and so on. 

In order to see he reason of the disagreement, we plot the form factors at a fixed momen- 
tum configuration ap^ = (0, 0, 0) and afc^r = (1, 0, 0) as a function of the light quark mass in 
Figure [Tl|. While we find a good agreement for f2{v ■ k^^), our results for /i(t> ■ /c^) + f2{v ■ k^,) 
is significantly lower than the Fermilab data. Furthermore, in the fit of the form ( |6.12|) the 
chiral limit of our data is lower than the data at finite arrig as shown in the plot, in con- 
trast to the Fermilab data, for which the chiral limit becomes even higher due to a positive 
curvature. 

We note that the renormalization of the vector current is made using the non-perturbative 
Z factors of heavy-heavy and light-light current in the Fermilab analysis A correction is 



then made perturbatively for the heavy-light current. Since our results are obtained with an 
entirely perturbative matching, systematic errors may enter differently. The effect of such 
a 'partial' non-perturbative renormalization for the NRQCD action is an issue of future 
investigation. 

We should also note that in Figure |Tl| the statistical error in our calculation seems much 
larger than that in the Fermilab data, despite much larger statistics in our calculation. We 
suspect that the main reason for the large statistical error in our data is a larger temporal 
extent of our lattice A't=48 compared to Nt=32 in the Fermilab work. The large temporal 
size and the large distance between t^^ and in our simulation renders the extraction of the 
ground state contribution very convincing as shown in Figure ^ which seems much better 
than the equivalent plot in |]TT|, but at the same time the statistical noise exponentially 



grows as the heavy-light meson evolves in the temporal direction |31,32 



B. and fiq^) 

A comparison of the form factors in the conventional definition f^{q^) and f^{q^) is 



made in Figure |T2|. Results from recent lattice calculations by APE 0, UKQCD [y] and 
Fermilab |]TT| are shown in the plot together with our data. 

We find that all data is consistent with each other for f^{q^), while our result is somewhat 
lower for Since /"(g^) is proportional to fi{v- k-„) + f2{v- /c^) up to a small correction 

of 0(f ■ /ETr/ms), the disagreement with the Fermilab result is the same one as we discussed 
in the previous subsection. 

The results of other two groups, APE and UKQCD, are lower than the Fermilab re- 
sult to but still higher than ours. We note that in their approach an extrapolation 
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in 1/Mp is necessary to predict the B meson form factor from the simulation results 
for lighter heavy quarks. Figure |13| shows such an extrapolation. The magnitude of 
$o(= {as{Mp)/as{MB)-^^^^f^/M^) in their results agrees with ours, but the APE re- 
sults show a negative slope in contrast to a fiat 1/M dependence of our data, leading to the 
APE value at the physical point considerably higher than ours. In the relativistic approach, 
the discretization error may be magnified toward heavier quarks, since the discretization 
error scales as a power of aM. Therefore, the dependence on the heavy quark mass can be 
badly distorted. Furthermore, the heavy quark expansion becomes questionable for lighter 
heavy quarks, and the extrapolation with a linear or quadratic function in 1/Mp may not be 
sufficient. Such effects are difficult to incorporate in the extrapolation, and the systematic 
error can be underestimated. 



VIII. PHENOMENOLOGICAL IMPLICATIONS 

A. Differential decay rate 

The differential decay rate of the semileptonic n^l^ui decay is proportional to the 

form factor f^{q^) squared. 



Therefore, if a reliable calculation of the form factor is available from lattice QCD, the ex- 
perimental data can be used to extract the CKM element \Vuh\- Our result for the differential 



decay rate divided by \Vub\ is listed in Table XIII and shown in Figure 14. 



The momentum configuration where data is available is limited to a large region, 
18 GeV^ < q^ < 21 GeV^, which corresponds to small recoil momenta. In the region 
above 21 GeV^ there is no data point because of large pion mass in the lattice calculations. 
However, the pole dominance near zero recoil region ( |2.11| ) and ( |2.12| ), which is confirmed 



in part by our lattice calculations, should become even better approximation in that region. 
Therefore, the theoretical uncertainty is under control in that large q^ region. A strategy 
to determine \Vub\ is, then, to measure the decay rate in the large q^ region, > 18 GeV^, 
and to use the lattice result 



F 



247r3 



j^^ = 1-18 ±0.37 ±0.08 ±0.31 psec~^ (8.2) 

JlSGeV^ 



The first error is statistical, the second is perturbative and the last error is the error from 
discretization and chiral extrapolation. 



B. L» ^ Tilv and D Klv 



As we found in Section |VI G| , the 1/M correction to the HQET form factors fi{v ■ k.^) 
and /2(f ■ /Ctt) is small. Although our data is only available for large heavy quark mass M > 
3.2 GeV and the charm quark mass is not covered, the result suggests that the semileptonic 
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decays of D mesons, D ttIu and D 
proposed by Burdman et al. 
The idea of 



Klv, may be used to constrain the form factors, as 



12 



12 



D 



is to consider a ratio of the differential decay rates of -B — vr/z/ and 
T[lv at a fixed recoil energy v ■ k.„, then the heavy quark symmetry tells us that the 
ratio is unity at the leading order, and the ratio of CKM elements |\4fe/V^d| may be extracted 
without model dependence. The method is, however, not so useful unless the size of 1/M 
(and higher order) corrections is reliably estimated. The lattice calculation could be used 
to evaluate them, as we attempt in this work. 

In the lattice calculation, the bulk of systematic errors, especially uncertainty in the 
perturbative renormalization, is canceled in the ratio of form factors with different heavy 
quark mass. This idea was extensively used by the Fermilab group ||33|j3^ in the lattice 
study of heavy-to- heavy decay, namely B D^*^lu, in which the heavy quark symmetry 
predicts a stronger constraint and the form factor is even normalized in the zero recoil limit 
up to a correction of 0(1/M^) that can be calculated on the lattice. The Fermilab group also 
considered the ratio for the heavy-to-light decay |]lT|. They calculated the form factors at b 
and c quark masses, and found a small but significant mass dependence in the HQET form 
factors, which might conflict with our findings. It is, therefore, important to extend our work 
toward lighter heavy quarks in order to investigate how the form factors are modified by 
the 1/M corrections. We also note that for this purpose the non-relativistic interpretation 
of the relativistic lattice action |jlOi employed in |[ll|] is best suited, because lighter heavy 



quark can be treated without large systematic errors. 



IX. CONCLUSIONS 

In this paper, we calculated the form factors and the differential decay rate for B ttIu 
on the quenched lattice using the NRQCD action. In the HQET form factors fi{v ■ A;,r) and 
f2{v ■ /c,r), the heavy quark mass dependence appears only in the form of the 1/M expansion. 
From calculations at several different heavy quark masses we found that the 1/M correction 
is not significant for these form factors. We found that the B* pole contribution dominates 
f^{q^) for small pion recoil energy. We also showed that the extrapolation to the soft pion 
limit suffers from large systematic errors, so that the discrepancies between P{q^ax) 
/b / fir in the soft pion relation, as seen, in the present simulation is not a serious problem. 

In order to avoid the model dependence, we did not assumed any particular functional 
form for the form factors. Instead, we carried out an interpolation in the region where our 
data are available. Although the accessible q^ region is rather limited, the prediction from 
chiral effective lagrangian may be used to extend the prediction toward q^axi ^'^^ ob- 
tained a partially integrated differential decay rate in the region 18 GeV< q^ < q^ax- We 

obtained — ^ / /^(g ) = 1.18 ± 0.37 ± 0.08 ± 0.31psec~^ where the first er- 

2AtC'^ JlSGeV^ 

ror is statistical error the second is the error from perturbative calculation and the third is 
the systematic error from the discretization and chiral extrapolation. 

The discretization error of O(a^) and the perturbative error are sizable. The first error 
can be reduced by performing simulation at several different lattice spacings and/or using 
different lattice actions. The reduction of the second error is more demanding. We need 
a non-perturbative renormalization to remove it. Another important source of uncertainty. 
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which wc did not include, is the quenched approximation, whose effect can be estimated 
only with simulations including dynamical quarks. We are planning future studies in these 
directions. 
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TABLES 



aMo 




n 


(0) (1) 

Pv! Pvi 




(2) 

Pvi 


10.0 




2 


-0.562 -0.572 




-0.421 


5.0 




2 


-0.554 -0.571 




0.205 


3.0 




2 


-0.540 -0.582 




0.446 


2.1 




3 


-0.529 -0.604 




0.559 


1.3 




3 


-0.509 -0.629 




0.657 


TABLE I. Renormalization constants for V4. 


aMo 


n 




(0) (1) (2) 

Pvl Pvl Pvl 


(3) 

Pvl 


(4) 

Pvl 


10.0 


2 




-1.250 0.366 13.705 


0.983 


1.047 


5.0 


2 




-1.087 0.232 4.678 


0.881 


0.977 


3.0 


2 




-0.915 0.091 1.605 


0.774 


0.893 


2.1 


3 




-0.772 -0.049 0.594 


0.690 


0.812 


1.3 


3 




-0.546 -0.235 0.188 


0.587 


0.668 


TABLE II. Renormalization constants for T4. 


aMo 






ah 




bh 


10.0 






0.16 




1.50 


5.0 






0.28 




1.12 


3.0 






0.29 




1.07 


2.1 






0.30 




1.06 


1.3 






0.31 




1.06 


TABLE III. Smearing parameters for the heavy-light meson. 
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id. Pb K q_ 

pOOO.qOOO (0,0,0) (0,0,0) (0,0,0) 

plOO.qlOO (1,0,0) (0,0,0) (1,0,0) 

plOO.qOOO (1,0,0) (1,0,0) (0,0,0) 

plOO.qllO (1,0,0) (0,-1,0) (1,1,0) 

pOOO.qlOO (0,0,0) (-1,0,0) (1,0,0) 

pl00.q200 (1,0,0) (-1,0,0) (2,0,0) 



TABLE IV. List of momenta for which the three-point correlator is measured. 
Three-momentum is given in the unit of 2Tr/La. 'id.' will be used in the tables of numerical 
results. 
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oMo=5.0 


id. vP 


vp 


vP 


vP 


vP 


vP 


vP 


vP 



pOOO.qOOO 
plOO.qlOO 
plOO.qOOO 
plOO.qllO 
pOOO.qlOO 
pl00.q200 



2.059(44 
2.052(52 
1.79(17) 
1.645(64 
1.554(56 
1.34(11) 



0.1874(53) -0.1874(53) 

0.2304(72) -0.1862(60) 0.010(10) 

-0.023(28) 0.023(28) 0.876(90) 

0.272(15) 0.0602(86) -0.740(30) 

0.216(13) 0.0603(97) -0.751(29) 



-0.660(17) 0.001(3) 

-0.496(53) 0.013(20) 

-0.121(26) 0.004(5) 

-0.063(23) 0.011(7) 



-0.713(18) -0.053(4) 
-0.241(26) 0.241(26) 
-0.358(20) -0.266(10) 
-0.321(19) -0.268(11) 



0.487(43) 0.073(20) -0.697(62) -0.547(70) 0.042(16) -0.765(71) -0.263(25) 



aMo=3.0 



'W 



id. 



V 



pOOO.qOOO 
plOO.qlOO 
plOO.qOOO 
plOO.qllO 
pOOO.qlOO 
pl00.q200 



2.087(33 

2.066(39 

1.81(16) 

1.621(57 

1.547(50 

1.343(94 



0.217(41) -0.2173(41) 
0.2802(60) -0.2126(47) 0.031(8) 
0.008(23) -0.008(23) 0.877(80) 
0.302(13) 0.0426(66) -0.704(26) 
0.231(11) 0.0390(73) -0.710(25) 
0.495(34) 0.073(15) -0.627(46) 



-0.632(12) 0.016(2) 

-0.416(44) 0.003(16) 

-0.152(22) 0.034(4) 

-0.096(18) 0.0393(53) 

-0.544(53) 0.076(13) 



-0.697(13) -0.081(3) 
-0.207(22) 0.207(22) 
-0.351(17) -0.274(9) 
-0.325(16) -0.268(9) 
-0.762(57) -0.296(20) 



aMo=2.1 



W 



7T 



w 



id. 



pOOO.qOOO 
plOO.qlOO 
plOO.qOOO 
plOO.qllO 
pOOO.qlOO 
pl00.q200 



2.108(30 
2.072(35 
1.85(16) 
1.607(51 
1.555(47 
1.321(86 



0.2552(41) -0.2552(41) 
0.3354(61) -0.2469(46) 0.057(7) 
0.050(22) -0.050(22) 0.898(73) 
0.341(12) 0.0172(58) -0.683(23) 
0.254(10) 0.0129(65) -0.686(22) 
0.504(32) 0.067(15) -0.556(41) 



-0.606(10) 0.030(2) -0.680(11) -0.105(3) 

-0.347(39) -0.008(14) -0.178(19) 0.178(19) 

-0.184(18) 0.060(4) -0.350(15) -0.278(8) 

-0.130(15) 0.064(5) -0.333(14) -0.268(8) 

-0.551(48) 0.117(16) -0.749(51) -0.315(19) 



aMo=1.3 



W 



IT) 



W 



TT 



id. 



V 



pOOO.qOOO 
plOO.qlOO 
plOO.qOOO 
plOO.qllO 
pOOO.qlOO 
pl00.q200 



2.140(25 
2.066(34 
1.94(16) 
1.519(45 
1.581(45 
1.278(71 



0.3278(40) -0.3278(40) 
0.4323(74) -0.3087(54) 0.103(6) 
0.136(24) -0.136(24) 0.943(66) 
0.399(14) -0.0227(63) -0.636(23) 
0.304(10) -0.0361(63) -0.655(19) 
0.551(30) 0.049(11) -0.478(31) 



-0.556(10) 0.057(2) -0.644(10) -0.146(3) 

-0.229(34) -0.030(12) -0.130(16) 0.130(16) 

-0.234(16) 0.103(5) -0.338(13) -0.269(8) 

-0.189(14) 0.110(5) -0.348(13) -0.269(7) 

-0.564(39) 0.180(14) -0.733(41) -0.349(17) 



TABLE V. Matrix elements at k=0. 13630. The first column represents the momentum 



configuration as defined in Table IV 



27 



aMo=5.0 


id. 

















pOOO.qOOO 
plOO.qlOO 
plOO.qOOO 
plOO.qllO 
pOOO.qlOO 
pl00.q200 



2.222(61) 

2.211(70) 

1.76(30) 

1.71(11) 

1.50(10) 

1.31(19) 



0.2754(88) 

0.326(11) 

0.051(40) 

0.338(31) 

0.262(26) 

0.556(78) 



-0.2754(88) 

-0.2756(99) 

-0.051(40) 

0.036(15) 

0.038(17) 

0.074(30) 



-0.005(13) -0.708(23) -0.009(4) 

1.11(15) -0.523(82) 0.071(28) 

-0.843(59) -0.126(46) -0.056(10) 

-0.826(56) -0.037(41) -0.037(13) 

-0.81(11) -0.53(12) -0.008(26) 



-0.772(24) -0.056(5) 
-0.224(41) 0.224(41) 

-0.408(36) -0.243(19) 
-0.322(34) -0.249(20) 
-0.77(12) -0.229(40) 



aMo=3.0 


id. vP 


vP 


vp 










y(^) 



pOOO.qOOO 

plOO.qlOO 
plOO.qOOO 
plOO.qllO 
pOOO.qlOO 
pl00.q200 



2.242(48) 

2.217(54) 

1.81(27) 

1.667(97) 

1.485(89) 

1.32(15) 



0.3050(72) 

0.375(10) 
0.063(31) 
0.362(27) 
0.268(21) 
0.564(59) 



-0.3050(72) 

-0.3003(80) 

-0.063(31) 

0.020(12) 

0.020(13) 

0.062(21) 



0.019(11) -0.679(17) 0.009(3) -0.755(18) -0.085(4) 

1.06(12) -0.408(65) 0.058(19) -0.175(32) 0.175(32) 

-0.789(50) -0.151(38) -0.020(7) -0.390(30) -0.256(15) 

-0.771(47) -0.076(32) -0.004(10) -0.326(29) -0.247(15) 

-0.726(81) -0.519(87) 0.034(18) -0.770(93) -0.285(30) 



ai\ /(,=2.1 


id. vr 


vP 


vP 


vP 


^(1) 


vP 







pOOO.qOOO 2.254(42) 0.3434(67) -0.3434(67) 

plOO.qlOO 2.214(50) 0.431(10) -0.3345(81) 0.050(10) 

plOO.qOOO 1.85(25) 0.092(30) -0.092(30) 1.05(11) 

plOO.qllO 1.634(86) 0.399(25) -0.005(10) -0.765(43) 

pOOO.qlOO 1.495(83) 0.287(20) -0.003(11) -0.743(42) 

pl00.q200 1.31(12) 0.590(51) 0.055(17) -0.679(61) 



-0.650(15) 0.024(2) 
-0.330(56) 0.046(16) 
-0.183(32) 0.013(6) 
-0.115(28) 0.024(8) 
-0.539(71) 0.076(15) 



-0.737(16) -0.111(4) 
-0.141(27) 0.141(27) 
-0.382(27) -0.263(13) 
-0.337(26) -0.247(13) 
-0.777(78) -0.314(26) 



aMo=1.3 


id. vr 


vP 


vP 






vP 


vP 





pOOO.qOOO 2.276(35) 0.4170(66) -0.4170(66) 

plOO.qlOO 2.193(43) 0.527(10) -0.3956(80) 0.1038(80) -0.595(12) 0.0531(20) 

plOO.qOOO 1.91(23) 0.146(35) -0.146(35) 1.02(10) -0.215(47) 0.019(17) 

plOO.qllO 1.486(72) 0.443(25) -0.041(11) -0.701(40) -0.225(27) 0.0644(74) 

pOOO.qlOO 1.529(76) 0.332(19) -0.049(11) -0.711(34) -0.182(23) 0.0754(82) 

pl00.q200 1.28(10) 0.620(47) 0.037(14) -0.575(45) -0.559(56) 0.149(15) 



-0.697(13) -0.1557(37) 
-0.098(23) 0.098(23) 

-0.352(23) -0.254(11) 
-0.356(22) -0.250(12) 
-0.765(62) -0.355(23) 



TABLE VI. Matrix elements vlf^ at «;=0. 13711. 
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aMo=5.0 



id. v}^^ 




vP 












pOOO.qOOO 2.433(92) 


0.382(15) 


-0.382(15) 












plOO.qlOO 2.42(10) 


0.448(21) 


-0.385(17) 


-0.031(21) 


-0.771(34) 


-0.020(60) 


-0.850(37) 


-0.059(8) 


plOO.qOOO 1.75(55) 


0.118(75) 


-0.118(75) 


1.47(31) 


-0.57(15) 


0.170(62) 


-0.197(77) 0.197(77) 


plOO.qllO 1.74(19) 


0.390(63) 


0.019(29) 


-0.94(11) 


-0.130(85) 


-0.120(23) 


-0.461(65) 


-0.203(35) 


pOOO.qlOO 1.33(17) 


0.292(50) 


0.014(32) 


-0.85(11) 


0.007(75) 


-0.081(28) 


-0.288(63) 


-0.224(36) 


pl00.q200 1.36(33) 


0.70(16) 


0.089(53) 


-1.06(24) 


-0.55(22) 


-0.061(53) 


-0.83(21) 


-0.216(74) 


aMo=3.0 


id. vr 
















pOOO.qOOO 2.439(71) 


0.411(12) 


-0.411(12) 












plOO.qlOO 2.413(78) 


0.494(17) 


-0.406(13) 


-0.001(17) 


-0.738(24) 


-0.000(4) 


-0.830(27) 


-0.091(6) 


plOO.qOOO 1.76(50) 


0.111(60) 


-0.111(60) 


1.38(25) 


-0.40(11) 


0.145(45) 


-0.127(60) 0.127(60) 


plOO.qllO 1.69(16) 


0.407(53) 


0.003(22) 


-0.864(96) 


-0.142(68) 


-0.076(16) 


-0.430(54) 


-0.230(28) 


pOOO.qlOO 1.31(15) 


0.283(42) 


0.004(24) 


-0.785(91) 


-0.043(58) 


-0.038(19) 


-0.294(53) 


-0.218(27) 


pl00.q200 1.37(26) 


0.70(12) 


0.057(38) 


-0.92(16) 


-0.55(15) 


-0.006(34) 


-0.83(16) 


-0.284(56) 


ai\ /(,=2.1 


id. 

















pOOO.qOOO 2.438(61) 0.449(11) -0.449(11) 

plOO.qlOO 2.392(72) 0.548(17) -0.439(13) 0.037(14) 

plOO.qOOO 1.81(46) 0.132(59) -0.132(59) 1.33(21) 

plOO.qllO 1.63(14) 0.441(49) -0.023(19) -0.837(83) 

pOOO.qlOO 1.33(14) 0.299(39) -0.014(22) -0.760(79) 

pl00.q200 1.35(22) 0.72(10) 0.047(31) -0.85(12) 



-0.705(22) 0.017(3) 
-0.305(97) 0.125(37) 
-0.169(57) -0.036(12) 
-0.088(49) -0.009(15) 
-0.572(12) 0.043(27) 



-0.807(24) -0.119(5) 
-0.086(51) 0.086(51) 
-0.410(47) -0.244(24) 
-0.310(47) -0.217(24) 
-0.85(13) -0.315(47) 



aMo=1.3 


id. 


vr 




vr 


vl'^ 


vr 


vl'^ 


vr 


vr 


pOOO.qOOO 


2.445(51) 


0.525(11) 


-0.525(11) 












plOO.qlOO 


2.350(61) 


0.643(17) 


-0.499(12) 


0.102(12) 


-0.643(17) 


0.050(3) 


-0.762(20) 


-0.169(5) 


plOO.qOOO 


1.86(41) 


0.163(68) 


-0.163(68) 


1.21(19) 


-0.190(83) 


0.088(35) 


-0.049(44) 0.049(44) 


plOO.qllO 


1.43(11) 


0.476(47) 


-0.058(21) 


-0.751(72) 


-0.200(46) 


0.024(12) 


-0.355(40) 


-0.242(20) 


pOOO.qlOO 


1.38(13) 


0.341(37) 


-0.054(22) 


-0.742(63) 


-0.160(41) 


0.042(14) 


-0.337(40) 


-0.222(21) 


pl00.q200 


1.35(17) 


0.741(86) 


0.023(26) 


-0.711(87) 


-0.610(96) 


0.126(25) 


-0.85(10) 


-0.367(42) 






TABLE VII. Matrix elements V^i'^ at k= 


=0.13769. 
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aMo 










5.0 


2.78(15) 




0.551(31) 


-0.551(31) 


3.0 


2.76(11) 




0.578(23) 


-0.578(23) 


2.1 


2.74(10) 




0.615(22) 


-0.615(22) 


1.3 


2.721(85) 




0.690(20) 


-0.690(20) 


TABLE VIII 


1\/T>^t,nv plpTTlPTlt, 1/^. 

±VXCtuX±./V C'lv/J.J.J.VjJ.l U ' /I 


^ at K=0. 13816 for the zero recoil configuration (pOOO.qOOO). 


n 


0.13630 


0.13711 


0.13769 


0.13816 


(0,0,0) 


0.48816(62) 


0.40756(68) 


0.34005(78) 


0.2723(11) 


(1,0,0) 


0.6209(28) 


0.5578(44) 


0.5065(73) 


0.459(15) 


(1,1,0) 


0.7424(43) 


0.6967(71) 


0.654(15) 


0.626(32) 


TABLE TX 


rinn pnpi^o"\7" n h! ( 1 


for spatial momenta (0,0,0), (1,0,0), (1,1,0) 


in unit of In /La. 


K 


0.13630 


0.13711 


0.13769 


0.13816 


aMo=5.0 


5.5702(21) 


5.5476(29) 


5.5318(39) 


5.5203(58) 


aMo=3.0 


3.6606(15) 


3.6376(20) 


3.6212(27) 


3.6088(40) 


aMo=2.1 


2.7992(12) 


2.7758(16) 


2.7588(22) 


2.7455(32) 


aMo=1.3 


2.0301(10) 


2.0058(13) 


1.9883(17) 


1.9741(25) 



TABLE X. Heavy-light meson mass uMb evaluated using ( |6.3| ) with av{l/a). 
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oMo=5.0 




0.13630 


0.13711 


0.13769 


0.13816 


(0,0,0) 


0.5887(21) 


0.5661(29) 


0.5503(39) 


0.5388(58) 


(1,0,0) 


0.6062(15) 


0.5821(22) 


0.5668(29) 


0.5554(40) 


aMo=3.0 




0.13630 


0.13711 


0.13769 


0.13816 


(0,0,0) 


0.5830(15) 


0.5600(20) 


0.5436(27) 


0.5312(40) 


(1,0,0) 


0.6071(13) 


0.5826(17) 


0.5665(23) 


0.5546(32) 


aMo=2.1 




0.13630 


0.13711 


0.13769 


0.13816 


(0,0,0) 


0.5750(12) 


0.5515(16) 


0.5346(22) 


0.5212(32) 


(1,0,0) 


0.6054(12) 


0.5804(15) 


0.5638(20) 


0.5512(28) 


aMo=1.3 




0.13630 


0.13711 


0.13769 


0.13816 


(0,0,0) 


0.5571(10) 


0.5328(13) 


0.5152(17) 


0.5010(25) 


(1,0,0) 


0.5984(12) 


0.5711(18) 


0.5532(24) 


0.5387(35) 



TABLE XI. Binding energy of the heavy-light meson a£'bin(p_B) for spatial momenta (0,0,0) 
and (1,0,0) in unit of 2-k jLa. 



^(000) 
'-^1+2 


^(010) 
•-"1+2 


^(100) 
•-^1+2 


^(001) 
•-"1+2 


^(011) 


^(002) 
•-"1+2 


0.413(74) 


3.79(97) 


-0.019(31) 


-0.53(66) 


-3.3(13.8) 


0.7(2.3) 


0.392(75) 


3.94(99) 


0.070(29) 


-0.59(66) 


-3.5(13.8) 


0.7(2.3) 


^(000) 


^^010) 


^^100) 


^^001) 






0.311(47) 
0.347(50) 


1.06(1.11) 
0.99(1.14) 


0.035(37) 
-0.020(37) 


-0.06(40) 
-0.04(40) 







TABLE XIL Global fit parameters in the form (|6.12| )-( |6.13| ). In each column, top and bottom 
numbers correspond to the result with ay(l/a) and ayijx jd). 
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V ' kj^ 


^2 

q 


Jl{V ■ Kn) + h[V ■ K7rJ 


hyv ■ f^w) 


J VI ) 


J [Q ) 


1 /It/ |2 j-p / J^2 




(GeV) 






(GeV ' ) 






(psec GcV j 




U.14>}5 


26.37 


A AO /OO \ 

0.98(23) 


A A f / 0A\ 

0.45(20) 


A /I / 1 \ 

0.84(18) 


'7 /I / 1 \ 

7.4(3.1) 


A AA1 7 / 1 r \ 

0.0017(15) 




u.iyio 


25.87 


0.95(20j 


0.44(19) 


n on/1 ci\ 
0.80(16) 


r r/o 0\ 


n noi /I '7\ 
0.021(17) 




(J. 2392 


25.36 


A A1 / 1 0\ 

0.91(18) 


0.44(18) 


0.76(15) 


4.4(1.6) 


A A/10/0"l\ 

0.042(31) 




(J.z87U 


24.86 


0.87(17) 


A A A / 1 '7\ 

0.44(17) 


A 70 { '\ A\ 

0.73(14) 


'7/ 1 o\ 

3.7(1.2) 


A A/^O / /I 1 \ 

0.062(41) 




n o o /I o 

(J.3o48 


24.35 


A 00/1 

0.83(17) 


0.44(15) 


A ^A { '\ A\ 

0.70(14) 


1 / AO \ 

3.18(98) 


A AO 1 / CT A\ 

0.081(50) 




L).oo27 


23.85 


0.79(17) 


A A A / 'I A\ 

(J. 44(14) 


A C'd { '\ A\ 

0.66(14) 


'70 / ^70 \ 

2.78(78) 


A AAA { TiCi\ 

0.099(56) 




0.4305 


OO O A 

23.34 


A T/^ / 1 0\ 

0.76(18) 


A /I / 1 0\ 

0.43(13) 


A ^0/1 A\ 

0.63(14) 


A£} i c: a \ 

2.46(64) 


0.115(60) 




C\ A 'TOO 


OO O A 

22.84 


A 'TO / 1 \ 

0.73(18) 


A /f / 1 0\ 

0.43(12) 


0.61(15) 


01 / c o\ 

2.21(53) 


A 1 01 ( 

0.131(62) 




U.OZDZ 


OO o o 

22.66 


A A / 1 \ 

0.69(18) 


A /I / 1 1 \ 

0.43(11) 


0.58(15) 


f\r\{ A A\ 

2.00(44) 


A 1 Ad { a A\ 

0.146(64) 




0.574U 


O 1 o o 

21.83 


0.66(18) 


A /I /" 1 A\ 

0.43(10) 


A r CT / 1 a\ 

0.55(14) 


1 00 /" 0'7\ 

1.82(37) 


0.161(65) 




O.D2io 


1 OO 

21.62 


0.64(18) 


A /100/AO\ 

0.428(93) 


A c / 1 /I ^ 
0.53(14) 


1.67(31 ) 


0.174(65) 




0.OD97 


OO OO 

20.82 


0.61(17) 


0.426(88) 


0.51(14) 


1.54(27) 


0.187(66) 


* 


0.7175 


on o 1 
20.31 


0.58(17) 


0.424(85) 


A AC\( '\ ^\ 

0.49(13) 


1 /I / /I \ 

1.43(24) 


A 1 A A f f*0\ 

0.199(68) 


* 


0.7d5o 


1 O O 1 

19.81 


0.56(15) 


A /100/'0/l\ 

0.422(84) 


A /I '7/ 1 0\ 

0.47(12) 


1 00/ oo\ 
\.i6\ll) 


A 01A/71\ 

0.210(71) 


* 


U.8132 


1 r\ OO 

19.30 


0.54(14) 


0.421(86) 


A /I CT / 1 1 \ 

0.45(11) 


1 /I /"o 1 \ 

1.24(21) 


A 00 1 / 7/^ \ 

0.221(76) 


* 


0.8d1U 


18.80 


A r / 1 \ 

0.52(13) 


A A C\/ A A\ 

0.419(90) 


A /iocr/Ao\ 

0.435(98) 


1 1 /oi \ 

1.16(21) 


A 001 / A \ 

0.231(84) 


* 


0.9088 


18.29 


0.50(12) 


0.417(97) 


0.421(92) 


1.09(21) 


0.240(94) 


* 


0.9567 


17.79 


0.48(11) 


0.42(10) 


0.407(92) 


1.03(22) 


0.25(11) 


* 


1.0045 


17.28 


0.47(13) 


0.41(11) 


0.40(10) 


0.97(23) 


0.26(13) 




1.0523 


16.78 


0.45(15) 


0.41(12) 


0.38(12) 


0.92(25) 


0.27(15) 




1.1002 


16.27 


0.44(19) 


0.41(14) 


0.38(15) 


0.87(27) 


0.27(17) 




1.1480 


15.77 


0.43(24) 


0.41(15) 


0.37(18) 


0.83(29) 


0.28(20) 





TABLE XIII. Numerical results for the form factors and for the differential decay rate. Error 

contains the statistical and estimated systematic uncertainties. Results with a * symbol in the 
last column are those obtained by interpolating the lattice data m. v ■ k-,^, while others involve an 
extrapolation. 
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FIGURES 



P3=(i,o,o), Mo.o.o) 




p,=(i,o,o),Mi.o.o) 








12 16 20 24 



12 16 20 24 



a! 



Pe=(1,0,0), k,=(0,-1,0) 



Pi 



Pe=(0,0,0),M-1,0,0) 



12 16 20 24 



12 16 20 24 



Pe=(1,0,0),M-1,0,0) 



^ 



12 16 20 24 



FIG. 1. Ratio i?^** (t; k-,r,PB for five combinations of fe^r and p^). Filled symbols represent the 
iio for vj^\ and open symbols are for Vi^\ Light quark is 
mass roughly corresponds to the b quark mass, i.e. aMo=3.0. 



ratio for vj^\ and open symbols are for Vi^\ Light quark is at «;=0. 13711, and the heavy quark 
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0.13711 
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0.13769 


A 


K = 


0.13816 



Q Q LJ I I I I I I I I I I I I I I I I I 

■ 0.0 0.1 0.2 0.3 

(ak f 

FIG. 2. Dispersion relation for pion. The lines represent the continuum form (|6.1| ). 




o 


K = 


0.13630 




□ 


K = 


0.13711 




o 


K = 


0.13769 




A 


K = 


0.13816 






0.2 



0.2 




o 


K = 


0.13630 




□ 


K = 


0.13711 




o 


K = 


0.13769 




A 


K = 


0.13816 





o 


K = 


0.13630 




□ 


K = 


0.13711 




o 


K = 


0.13769 




A 


K = 


0.13816 






(apj 



FIG. 3. Dispersion relation for the heavy-light meson at aMg = 5.0, 3.0, 2.1 and 1.3. The 
lines represent the non-relativistic form ( |6.2[) with perturbatively calculated meson mass aMs- 
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1.0 
0.8 



0.4 
0.2 



1 , 1 , 
'_ a'" [f/vkj+fjvkj ] 


1 

1 




? # J 


1,1, 


1 







0.2 



0.4 
<3 v-k 



0.6 



0.8 



71 



FIG. 4. A typical plot of the form factors fi{v ■ /c^) + f2{v ■ A;^) (open symbols) and f2{v ■ k^^) 
(filled symbols) in the lattice unit. Parameters are aMQ=3.0, k=0. 13630, and av(l/a) (circles) or 
ay(7r/a) (squares) are used for the perturbative matching. The data for av{TT/a) is slightly shifted 
in horizontal direction for clarity. 




a v-k 

n 

FIG. 5. The renormalization group invariant form factors <I>i+2(w ■ fcvr) (open symbols) and 
^2iv ■ Aitt) (filled symbols) for different values of uMq with fixed light quark mass k=0. 13630. 
Symbols denote the data at aAfo=5.0 (circles), 3.0 (squares), 2.1 (diamonds) and 1.3 (triangles). 
Solid and dashed lines show the fit ( |6.12 )-( |6.13|) for the heaviest (aMo=5.0) and the lightest 
(aMo=1.3) heavy quark masses respectively. 
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FIG. 6. Light quark mass dependence of fi{v ■ k^^) + f2{v ■ k^^) (open symbols) and f2{v ■ k-,^) 
(filled symbols) for aMo=2.1. Symbols denote the data at k=0. 13630 (circles), 0.13711 (squares), 
0.13769 (diamonds) and 0.13816 (triangles). Three thin solid lines, from above to below, show the 
fit ( |6.13| ) for f2(v ■ k-j^) with three k values, from heaviest to lightest. Four thin dashed lines, on 
the other hand, correspond to the fit (|6.12| ) for the data of fi{v-kT^) + f2{v ■ k^^) for four values of 
K,. Thick lines represent the limit of physical light quark mass. 
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FIG. 7. Form factors fi{v ■ /c^) + f2{v ■ k^^) (open symbols) and f2{v ■ k^^) (filled symbols) at 
physical mass parameters. The points with symbols are obtained by interpolation in • A;^, while 
others involve extrapolations. 
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FIG. 8. Soft pion limit of fo{q. 



2 1 

max J 



/■iriB 



[fi{v ■ + f2{v ■ K)] at aMo=3.0. The dashed 



line is a linear fit in (arrij^)'^, while the solid curve includes the term (am^). A result of the fit 
( 6.12| ) is given by an open triangle, which should be equal to /b/Ztt (filled triangle) in the soft pion 
theorem. 
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FIG. 9. Same as Figure but with estimated systematic errors. 
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FIG. 10. Form factors fi{v ■ k^^) + f2{v ■ k^^) (open symbols) and f2{v ■ k^^) (filled symbols) 
at physical mass parameters. Squares represent the results of while our data presented in 
Figure |^ is now plotted with gray symbols. 
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FIG. 11. Form factors fi{v ■ kj^) + f2{v ■ /c^) (open symbols) and f2{v ■ k^^) (filled symbols) for 
a fixed momentum configuration ap^ = (0,0,0) and afcjr = (1,0,0) are plotted as a function of 
light quark mass arUq/uQ. Triangles are results of the Fermilab group |]l^] for a heavy quark mass 
close to the h quark mass. Our results are shown for aM^ = 3.0 (circles), 2.1 (squares) and 1.3 
(diamonds). Squares and diamonds are shifted in the horizontal direction for clarity. Lines show 
the global fit (|1|) and (|6ll ). 
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FIG. 12. Comparison of the results for the form factors f^{q^) (filled symbols) and f^{q^) 
(open symbols). Data are from APE 0] (up trianlgles), UKQCD Q (down triangles) and Fermilab 
IjlTI (squares). Our results are plotted by circles and error bands are shown by dashed lines. 



1.5 - 



1.0 - 



0.5 - 



0. 



\7 



o o o o 



<> 



<> 



<> 



<> (GeV^^^) 



(GeV"^^^) 



0.2 0.4 0.6 
1/Mp (GeV^) 



0.8 



1.0 



FIG. 13. 1/Mp dependence of the form factors = {as{Mp)/as{MB)~^^^^ f+ /VW (filled 
symbols) $o = (a^(Mp)/a^(MB)-2/ii/°\/Mp f+iq"^) (open symbols) at a fixed v (= 0.845 
GeV). Simulation results from the APE collaboration are shown by diamonds, and their linear 
and quadratic extrapolation to the B meson mass is plotted by down and up triangles, respectively. 
Our results are given by circles. 
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q (GeV^) 

FIG. 14. Differential decay rate of the semileptonic B ttIu decay. The points with symbols 
are obtained by interpolation in v ■ k^j^, while others involve extrapolations. 
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